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SUMMARY 


The  problem  of  assessing  the  errors  in  the  parameters  obtained  from 
a  curve-fitting  process  is  considered,  and  a  scheme  which  may  be  applied 
toward  the  solution  of  such  problems  is  obtained.  ,This  method  is  then 
specialized  to  the  problem  of  finding  the  errors  in  the  calculated 
stability  parameters  of  an  airplane,  and  an  example  is  given. 


INTRODUCTION 


Curve-fitting  procedures  have  found  places  in  nearly  all  branches 
of  engineering;  in  particular,  the  aeronautical  engineer  may  apply  these 
methods  to  the  calculation,  from  flight  data,  of  the  stability  parameters 
of  an  airplane  (references  1  and  2).  Whether  least  squares  or  any  of 
the  profusion  of  graphical  methods  which  exists  is  used  for  this  curve- 
fitting  process,  questions  of  the  errors  in  the  calculated  parameters 
are  bound  to  arise . 

Although  there  is  a  considerable  amount  of  literature  on  the  subject 
of  least  squares  and  curve  fitting,  comparatively  little  is  to  be  found 
on  the  related  subject  of  errors.  What  literature  does  exist  (e.g., 
reference  3)  attacks  the  problem  from  the  point  of  view  of  statistics, 
arriving,  finally,  at  a  quantity  called  the  variance.  This  quantity, 
while  giving  a  satisfactory  reply  to  the  error  question  when  applied  to 
fitting  a  set  of  data  to  a  straight  line  when  only  one  measurement  is 
subject  to  error,  is  far  from  adequate  for  other  curve-fitting  problems. 
One  does  not  have  to  look  far  to  find  the  reason  for  this;  it  is  that  no 
method  has  as  yet  been  devised  for  calculating  the  variance  when  either 
the  fitted  curve  is  not  linear  or  when  more  than  one  measured  quantity 
is  subject  to  error.  This  latter  objection  is  not  pertinent,  perhaps, 
when  the  problem  of  calculation  of  stability  parameters  is  considered, 
for,  although  both  the  input  (control- surface  deflection)  and  the  response 
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may  be  subject  to  error,  it  is  frequently  assumed  that  only  the  output 
is  fallible.  This  approximation  is  particularly  good  when  only  free- 
oscillation  data  are  analyzed.  The  first  objection  is,  however,  more 
serious,  for,  assuming  even  that  all  quantities  remain  within  the  so- 
called  linear  range,  so  that  the  response  satisfies  a  linear  differential 
equation,  the  response  will  certainly  not  be  a  linear  function  of  the 
parameters . 

It  is  shown  in  the  body  of  this  report  that  the  solution  to  the 
error  problem  depends  upon  a  more  or  less  arbitrary  definition,  the 
only  criterion  as  to  the  choice  of  the  definition  being  the  usefulness 
of  the  solution  to  which  it  leads.  In  this  respect,  it  may  be  said  of 
the  method  derived  herein  that  it  is  not  difficult  to  apply  and  appears 
to  lead  to  reasonable  values  of  the  errors.  A  relationship  between  the 
error  formula  given  herein  and  the  classical  formula  for  the  variance 
(in  the  linear  case  when  such  a  formula  exists)  is  established. 


ANALYSIS 

Relation  of  the  Problem  to  Aerodynamics 


Suppose  one  has  a  set  of  data  which  represents  a  time  history  of,  .  ( 
for  example,  pitching  velocity  of  an  airplane  in  response  to  an  elevator 
deflection  S(t).  If  q(t)  represents  the  pitching  velocity,  then, 
under  certain  simplifying  assumptions,  it  is  shown  in  reference  1  that 
q(t)  satisfies  the  differential  equation 

£!a+-b  kq  =  cx  —  +  c06  (1) 

dt2  dt  <it 

where  the  constants  b,  k,  C-j_  and  C0  are  functions  of  the  stability 
derivatives  of  the  airplane. 

It  may  be  verified  by  differentiation  that  the  following  function 
represents  the  general  solution  of  equation  (l): 


-  o 


o( t) dv  e 


where  A^  and  A  2  are  constants 
(dq/dt)t=0,  and  5(o),  while  Xi 
istic  equation 


depending  on  the  initial  conditions  q(o), 
and  \z  are  the  roots  of  the  character - 


\2  +  bX  +  k  =  0 
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It  is  often  desirable  to  be  able  to  ascertain  the  "best"  values  of  the 
constants  b  and  k  (or  X^  and  Xg),  Co  and  Ci  corresponding  to  the  given 
time  histories  of  q(t)  and  S(t).  This  problem  is  the  subject  of  refer¬ 
ences  1  and  2. 

Suppose  then,  that  by  some  means  these  constants  have  been  evaluated. 
There  will  be  a  certain  "error"  in  these  values,  however,  due  to  the 
experimental  error  in  the  data,  the  simplifying  assumptions  mentioned 
above,  and  other  causes.  The  question  of  the  magnitudes  of  these  errors, 
given  a  certain  error  in  the  data,  is  considered  in  the  present  report. 


General  Discussion  of  the  Problem 


In  order  to  state  adequately  the  problem  of  errors,  it  is  first 
necessary  to  give  a  precise  statement  of  the  curve-fitting  problem.  To 
this  end,  consider  a  physical  quantity  qe(t)  which  is  measured  at  t  =  t0, 
tf,  tN,  where  tG<ti  <...<%.  It  is  assumed  that  various  theo¬ 

retical  considerations  would  indicate  that  qe(t)  should  be  one  of  the 
functions  of  the  set  q^x^Xs, . . .  ,xm)  ,•  that  is,  there  should  exist 
values  of  the  parameters  x-l,  x2,  ...,  xm  such  that 

Te(^i)  =  <l(‘ti»x1,x2,  ...,xm)  (3) 

for  all  i  =  0,  1,  . ..,  N.  However,  because  of  certain  unknown  errors 
in  q6,  equation  (3)  is  not  exactly  satisfied  for  all  i.  It  is  desira¬ 
ble  then  to  find  those  values  of  the  parameters  x^  which  "most  nearly" 
cause  equation  (3)  "to  be  satisfied.  One  means  of  doing  this  is  to  define 
the  "best"  values  of  x-l,  ...,  xm  as  those  values  of  the  parameters 
which  make 

N 

M  =  ^  [ Q.( , x^ , Xg f . . . , Xjq)  -  qe(tp)}  (^0 

i=o  . 

a  minimum.  The  process  of  minimizing  M  is  called  curve  fitting  by 
least  squares,  and  this  general  problem  is  considered  in  reference  2. 

As  for  the  error  problem,  a  careful  study  of  the  extant  literature 
on  statistics  will  lead  one  to  the  conclusion  that  the  "true"  value  of  a 
quantity,  upon  which  the  intuitive  definition  of  error  rests,  has  never 
been  meaningfully  defined.  Clearly,  however,  the  values  obtained  for  the 
errors  in  a  particular  problem  will  depend  upon  this  definition.  When 
the  usual  probabilistic  statement  of  the  error  problem  is  chosen,  the 
quantity  called  the  variance  arises,  along  with  the  objections  raised  in 
the  introduction  which  are  concomitant  with  it. 
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Various  nonprobabilistic  statements  of  the  problem  are  conceivable, 
and  the  statement  given  below  has  been  chosen  from  among  them  as  a  useful 
one.  These  nonprobabilistic  statements  are  based  on  the  concept  of  the 
"sensitivity"  of  the  function  (  q(t,x1, . . . ,xm)  with  respect  to  the  parame¬ 
ters.  That  is,  the  question  is  asked:  If  one  of  the  parameters  is 
changed  slightly,  does  q  also  change  only  slightly,  or  does  q  change 
by  a  large  amount?  If  a  small  change  in  a  parameter  ends  in  a  large 
change  in  q,  q  is  said  to  be  sensitive  to  changes  in  that  parameter, 
and  it  appears  clear  that  any  meaningful  theory  of  errors  should  lead  to 
relatively  small  errors  in  such  a  parameter.  On  the  other  hand,  if  q 
changes  by  only  a  small  amount  when  a  parameter  is  changed,  q  is  insen¬ 
sitive  with  respect  to  that  parameter,  and  the  theory  should  result  in  a 
large  error.  A  quantitative  discussion  follows. 


Statement  and  Solution  of  the  Problem 


In  place  of  the  elusive  term  "experimental  error,"  we  shall  intro¬ 
duce  the  concept  of  "residual,"  defined  by  the  equation 

€i  =  q(t^,x1,x2, . . . ,xm)  -  qg(ti)  (5) 


so  that  equation  (4)  may  be  written  in  the  equivalent  form 

N 


i=o 


Suppose  the  curve-fitting  problem  has  been  solved;  that  is,  suppose 

values  x^0),  x2(°),  ...,  xm(°)  of  the  parameters  x1,  x£,  ... ,  xm, 
respectively,  have  been  found  which  minimize  M.  Let 


(o) 


7  *2  '  9 


y 


(le^i) 


so  that  the  minimum  value  of  M  is 


N 


i=o 


One  possible  statement  of  the  error  problem,  which,  with  certain 
modifications,  will  be  used  in  this  report,  is  the  following.  Choose 
any  set  of  numbers  Axx,  Ax2,  . . • ,  Axm,  and  let 


xk(0)  +  ^k 


( k — 1,2,  •  •  •  ,  ni) 
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The  question  then  is:  How  large  may  the  quantities'  |Axk|  be  so  that 

t1,x1(1),x2(1),...,xm(1)  -  q  I  <|e±  <°)  |  (6) 


for  all  i  =  0,  1,  . . . ,  N?  This  maximum  value  of  |Axk|  may  then  be 
defined  as  the  allowable  error  in  xk.  This  problem  can  be  solved,  but 
various  objections  can  be  raised  as  to  its  significance  as  a  statement 
of  the  error  problem.  In  particular,  it  is  easy  to  see  that  this  estimate 
of  the  errors  may  be  far  too  optimistic,  for  if  only  m  out  of  the  total 
of  (N+l)  values  of  epf0)  are  zero,  then  the  only  way  in  which  inequal¬ 
ity  (6)'  may  be  satisfied  for  all  i  is  for  x^1)  =  xk(°)  for  all  k. 

That  is,  since  the  number  m  of  parameters  is  usually  far  less  than  the 
number,  (N+l),  of  data  points,  we  can  say  that  even  though  most  of  the 
residuals  be  arbitrarily  large,  the  result  is  that  the  error  in  the 
parameters  is  zero I 

This  objection  may  be  overcome  by  relaxing  the  requirement  (6), 
substituting  in  its  place  a  mean-square  inequality  obtained  by  squaring 
both  sides  of  (6)  and  summing  over  i.  The  resulting  inequality  will 
then  be  required  to  hold  instead  of  the  inequality  (6).  Symbolically, 
it  will  be  required  that 


Thus,  condition  (6)  is  made  to  hold  only  in  a  mean-squdre  sense  over  the 
whole  range  of  t,  while  it  may  not  be  true  for  some  particular  values 
of  tq.  An  allowable  error  in  xk  is  then  defined  as  any  value  of  Axk 
for  which  inequality  (7)  holds. 


This  last  could  reasonably  be  used  as  the  definition  of  the  error 
in  xk.  In  the  interest  of  ease  in  calculation,  however,  a  Taylor's 
series  expansion,  with  only  the  first-order  terms  being  retained,  will 
be  used, to  linearize  the  problem.  This  may  be  done  since  neither  the 
residuals  eq(°)  nor  the  errors  Axk  may  be  too  large;  if  they  were,  it 
may  be  said  that  there  is  something  wrong  with  the  theory  which  predicted 
that  equation  (3)  will  be  approximately  true,  or  with  the  experiment 
leading  to  the  data  qe(tp),  or  that  the  experiment  has  not  been  properly 
designed.  In  order  to  shorten  the  formulas,  the  following  notation  will 
be  used:1  " 

1For  the  sake  of  clarity,  it  should  be  mentioned  that  two  different  types 
of  subscript  are  used  in  this  report,  attached  to  the  same  quantity  q. 
Of  the  experimental  data  qe,  we  have  already  spoken;  in  addition,  the 
quantities  qk  are  now  defined  by  equations  ( 8) .  There  need  be  no 
confusion,  however,  for  the  subscript  e  will  always  be  used  to  denote 
the  data;  other  subscripts  will  be  used  as  defined. 
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^i)  —9.  )>*•*.>  xm^  ^ 


,,  x  dqCtijX-^0),...^0)]  _  rt  N 

9.k(^i)  —  r  >  ( k-l>2, . . . ,m) 

dxk 

We  shall  also  write  Aq(tp)  in  place  of  the  difference 


9  t ^ ^  j •  jXj^(^)  -  q  tp  ,  x^  • . . , Xj)j  ^ 

_  J  L  _ 

so  that  the  inequality  (7)  defining  the  errors  in  the  xk  becomes 

N  2 

^  Aq(t±)J  <M^o)  (7) 

i=o 

If  the  errors  in  the  parameters  are  not  too  large,  it  is  clear  that 
the  following  equation  is  approximately  true: 

m 

Aq(ti)  =  qk(ti)Axk 


Utilizing  this  linearization  process,  inequality  (7)  becomes 

R  T  m  Is 


^  ^  9k(ti)Axk  <  M(o) 


i=o  [_  k=i  J 

The  modifications  of  the  definition  of  errors  have  finally  been  completed. 
We  now  define  the  error  in  xk  as  the  largest  value  of  |Axk|  for  which 
inequality  (9)  holds,  regardless  of  the  values , of  the  other  Ax's. 

In  order  to  be  able  to  draw  some  pictures,  it  will  now  be  assumed 
that  m  =  2.  The  generalization  to  larger  m  will  be  presented  after¬ 
wards.  Thus,  it  is  assumed  that  the  function  q  is  dependent  on  two 
parameters  only: 

q  =  q(ti,x1,x2) 
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For  brevity  in  the  formulas,  we  shall  write  for  Ax^  and  shall 

define 

A 
N 


Hu  =  ^  [^1(^1)]“ 


1=0 

N 


<ll2  =  ^ 
i=o 
N 

0-22  =  ^  [l2(ti)]2 


1=0 


(10) 


The  inequality  (9)  defines  a  certain  domain  in  the  (gx,g2)  plane  whose 
boundary  is  given  by  the  following  equation: 

Hilda.2  +  2q12^1|2  +  Q.22^2  =  (ll) 


The  graph  of  equation  (ll)  is  either  an  ellipse,  a  parabola,  or  an  hyper¬ 
bola.  It  is  clear  that  if  the  definition  of  error  which  has  been  chosen 
is  significant,  equation  (ll)  must  represent  an  ellipse,  for  if  it  did 
not,  infinite  errors  would  be  obtained.  It  is  just  as  well  that  the 
proof  that  equation  (ll)  does  indeed  represent  an  ellipse  be  given  here, 
for  it  is  very  simple.  Consider  the  discriminant 

^11^22  "  ^12 2 

of  equation  (ll),  which  becomes,  using  equations  (lO), 


^11  ^22 


N  N  ■  r  N 

q122=  Y,  Mu))2  Y,  Z  Sidth/ti) 


1=0 


x=o 


L  1=0 


The  curve  in  question  is  an  ellipse  if  and  only  if  the  discriminant  is 
positive.  However,  the  inequality 


N 


N 


r  N 


7  [li(ti)  ]  7  Icl2('ti)3  " 


1=0 


1=0 


Hi(  ^i)  l2(^i) 


L  i=o 


>  0 
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is  precisely  the  ■well-known  Schwarz  inequality.  As  is  also  well  known, 
the  equality  cannot  occur  unless  qx(ti)  is  proportional  to  q£(ti)  for 
all  i,  and  as  this  implies  that  the  parameters  xx  and  x2  are  not 
independent,  it  can  he  assumed,  without  loss  of  generality,  that  q^^i) 
and  q.£(ti)  are  not  proportional.  This  completes  the  proof. 

If  the  ellipse  which  represents  equation  (ll)  should  he  drawn,  a 
graph  similar  to  the  one  shown  in  sketch  (a)  would  he  obtained. 

£2 


Since  there  are  no  first-order  terms  in  equation  (ll),  the  ellipse  must 
he  symmetric  with  respect  to  the  origin.  Interpreting  this  curve,  it 
may  he  said  that  any  point  which  lies  inside  the  ellipse  has  coordi¬ 
nates  which  define  an  allowable  error  in  x^  and  x2.  However,  in  the 
definition  of  the  error  in  the  parameters,  it  was  said  that  the  error 
in  xjj  is  the  largest  value  oflAx^l  for  which  inequality  (9)  holds, 
regardless  of  the  values  of  the  other  Ax’s.  For  this  reason,  in  order 


NACA  TN  2820 


9 


to  define  the  error  in  the  parameters  hy  means  of  the  diagram,  we  must 
enclose  the  ellipse  of  sketch  (a)  in  a  rectangle,  as  in  sketch  (h) . 


^2 


Any  point  lying  within  this  rectangle  is  considered  as  defining  an 
allowable  error.  If  Hi  and  H2  denote  the  errors  (i.e.,  the  maximum 
allowable  errors)  in  xx  and  x2,  respectively,  then  ,an(i  H2  are 
obtained  as  the  maximum  values  of  g1  and  lying  in  the  rectangle. 


If  m  =  2  as  in  the  above  example,  the  errors  can  always  he  found 
graphically.  However,  if  m>2,  an  analytical  method  must  he  used. 

Such  a  method  will  now  he  found  for  the  case  where  m  =  2,  and  the 
generalization  will  then  be  shown. 

It  is  clear  that  the  sides  of  the  rectangle  drawn  as  solid  lines 
in  sketch  (h)  are  the  tangents  to  the  ellipse  at  the  points  where 


Ik 

d*2 


Calculating  these  derivatives  from 


equation  (ll), 


we  obtain 
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and  the  quantities  qj^),  qj^tjJ  are  defined  in  equation  (8).  If  Eh 
denotes  the  maximum,  value  of  satisfying  equation  (l4),  so  that  Wh 

is  the  error  in  x^,  the  problem  is  to  find  Eh.  It  may  be  solved  by 
setting  the  (m-l)  derivatives  of  the  form 


( 2  , « .  < ,h-l,h+l , . . . ,m) 


equal  to  zero  and  solving  the  resulting  equations  along  with  equation  (l4) 
for  ^h*  0ne  need  not  go  through  this  process  each  time  the  errors  in 
a  problem  are  to  be  calculated,  for  it  can  be  done  once  and  for  all  as 
follows:  Differentiating  equation  (l4),  it  may  be  seen  that 


Sgh 


m 

I 


^  qhklk 

k=i 


Setting  these  derivatives  equal  to  zero,  we  obtain  the  following  (m-l) 
equations  which  are  to  be  solved  simultaneously  along  with  equation  (l4) 
for  !h: 


m 

y  ‘  C ,  j — 1,2 , . . .  ,h  — 1  ,h+l ,  •  «  • , m 


(16) 


k=x 


Equation  (l4)  can  be  written  in  the  form 

m  /  m  \ 

X  (  =  «(o)  (H) 

<5=1  \k=i  / 

Taking  equation  (l6)  into  account,  it  may  be  seen  that  the  parenthesized 
quantity  is  zero  except  when  j=h,  and  so  we  obtain  the  equation 


m 

I 

k=i 


fyik^k  = 


M 


(°)_ 

£h 


(17) 


This  equation  and  equation  (l6)  are  m  linear  equations  to  be  solved 
for  1^.  This  may  be  done,  in  general,  by  means  of  Cramer’s  rule  (the 
method  of  determinants). 


12 
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Suppose  that  the  given  approximating  function  q(t)  were  linear  in 
the  parameters,  so  that 

<l(ti)  =  qi1(t1)x1  +  q2(ti)x2  +  ...  +  q.m(ti)xm 

where  the  f -unctions  qj  (t-jJ  are  independent  of  the  parameters  xk. 
Differentiating  this  last  equation,  we  have  that 


dq.(tj) 


=  QkC^i)  (k=l,  ...  ,  m) 


and,  utilizing  equation  (15),  that 

N 

<ljk  “  Z 

i=o 

In  reference  3,  the  familiar  normal  equations  obtained  from  the  least- 
squares  process  are  shown  to  he 


where' 


m 

I 


N 


=  le,k 


(21) 


CL€, k  =  qe(ti)qk(ti) 


i=o 

and  qe(ti)  represents  the  given  data.  Let  D  denote  the  determinant 
of  the  coefficients  in  equation  (21),  and  let  he  the  minor  of  q^ 

in  D,  so  that  D  and  are  given  by  equations  (l8)  and  (19),  respec¬ 
tively.  Then,  if  a2  denotes  the  mean-square  error  in  the  data,  while 
0k2  denotes  the  mean-square  error  to  he  expected  in  the  parameter  xk 
(i.e.,  a2  and  ak2  are  the  variances  of  the  data  and  of  x^,  respec¬ 
tively),  it  is  shown  in  reference  3  that  if  q(t^)  is  a  linear  function 
of  the  parameters. 


ah 


V  D 


(22) 


"The  comma  is  placed  after  the  subscript  e 


different  sort  of  symbol  from  the  other  subscripts.  Thus,  j  and  k 


to  indicate  that  it  is 
Thus 

to  vary,  assuming  the  values  1,  ...,  m;  however,  the  subscript  e 
used  merely  to  indicate  that  the  experimental  data  qe(t^) 
the  definition  of  q 


are 
is 

is  used  in 


Le,k' 
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Noting  now  that  a2  affords  an  estimate  for  the  error  in  the  data  when 
the  statistical  conception  of  the  error  problem  is  considered,  while 

occupies  a  similar  position  in  the  present  theory,  a  relationship 
between  the  variance  in  a  parameter  and  what  has  herein  been  called  the 
allowable  error  in  that  parameter  can  be  established.  Comparing  equa¬ 
tions  (20)  and  (22),  it  may  be  seen  that  the  ratio  of  the  variance  in  a 
parameter  to  the  statistical  measure  of  the  error  in  the  data  is  the 
same  as  the  ratio  of  the  error  in  a  parameter  to  the  nonstatistical 
measure  of  the  error  in  the  data.  In  symbols, 


_  2 
~h 


(23) 


This  is  a  remarkable  fact  that  in  formulating  a  theory  of  errors  which 
abstains  from  the  use  of  probability  theory,  a  notion  of  error  has  still 
been  defined  which  bears  as  intimate  a  relation  to  the  statistical  ideas 
as  that  described  by  equation  (23).  Of  course,  this  relationship  has 
been  proved  only  in  .the  case  of  linear  curve-fitting  problems.  From, 
this,  however,  it  might  appear  reasonable  that  even  for  general  nonlinear 
problems  equations  (20)  and  (23)  (and,  therefore,  (22))  afford  an  expres¬ 
sion  for  the  variance.  This  leaves  an  interesting  problem  for  future 
research. 


There  is  an  important  conclusion  to  be  drawn  from  equation  (23) • 
If,  as  before,  N  +  1  is  the  number  of  data  points,  then  cr  and  M 
are  related  (reference  3)  by  "the  equation 


N+l 


Thus,  from  equation  (23), 

(N+l)ah2 

which  implies  that  is  larger  than  ah.  This  result  should  not  be 

too  surprising,  for  it  will  be  recalled  that  Oh2  is  a  measure  of  the 
probable  error  in  Xh,  while  we  have  tried,  in  *h>  to  define  some  sor 
of  maximum  error. 


APPLICATION  TO  AN  EXAMPLE 


Referring  to  equation  (2),  if  there  is  some  value  ta  of  time  such 
that  5(t)  =  0  for  all  t >ta  (i.e. ,  if  a  pulse  elevator  input  has  been 
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applied  to  the  airplane),  it  may  he  seen  that  for  t>ta,  the  pitching 
velocity  is  a  sum  of  two  exponentials 

q  =  +  B2el2t 


where  B1,B2,A1  and  X2  are  constants  depending  on  the 
stability  derivatives  of  the  airplane. 

In  general,  such  data  are  oscillatory;  that 
is,  11,12,Bi  and  B2  are  complex.  Let  X,x  =  Z  +  z’i, 

Bi=^  (p+3'i),  i2  =  -  1.  Since  q  is  real,  \s=l  - Z * i, 
B2  *  \  (P-P’l),  and 

q  =  e^(p  cos  Z't  -  p'  sin  Z't)  (24) 


The  data  given  on  this  page  are  actual  flight  data 
which  represent  the  pitching  velocity  of  a  test  airplane 
in  response  to  an  elevator  deflection  which  was  zero  for 
all  t >0.4  (t  =  0  is  taken  at  the  beginning  of  the  pulse). 
These  data  are  fitted  to  a  sum  of  exponentials  in  refer¬ 
ence  2,  where  the  results 

« 

z  =  -1.366  p  =  o.6i4 

=  3*071  p'  =  -0.208  ,  (25) 

m(°)  =  0.000895 

.  v 

were  obtained.  The  errors  in  these  parameters  will  now  be 
calculated,  utilizing  equation  (20)  and  substituting  z,  Z*,  P, 
and  p'  for  x1,x2,x3,  and  x4,  respectively. 


■ti 

qe(ti) 

0.4 

0.224 

.5 

.120 

.6 

.020 

*7 

-.057 

.8 

-.112 

*9 

-.148 

1.0 

-.160 

1.1 

-.150 

1.2 

-.127 

1*3 

-.097 

1.4 

-.062 

1*5 

-.032 

1.6 

-.005 

1*7 

.017 

1.8 

•030 

1.9 

.036 

2.0 

.035 

2.1 

.032 

2.2 

.027 

2.3 

.020 

2.4 

.015 

2.5 

.011 

2.6 

•  008 

2.7 

.005 

2.8 

.003 

2.9 

.001 

3.0 

0 

3.1 

0 

3.2 

0 

Calculation  of  the  errors  in  l,  i\  3,  andP'.-  Taking  derivatives 
from  equation  (24),  one  may  see  that 


8q 

8Z 

=  teZt(p 

cos  Z't  - 

3’ 

sin  Z't) 

8q 

bV 

=  -te  ^(p 

sin  Z’t  + 

3' 

cos  Z't) 

8q 

=  e  ^  cos 

Z't 

8p 

8q 

=  -e^t  sj_n 

Z't 

dp* 

(26) 
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Identifying  l  with  the  parameter  xlf  V  with  x2,  (B  with  x3, 
and  with  x4,  it  may  he  seen  that  q.1(ti)  £=  — j  is  given  in 

column  11  of  table  I  and  q2(ti)  is  given  by  minus  column  15,  while 
‘Igf'ti)  and  q.4(ti)  are  given,  respectively,  by  column  6  and  by  minus 
column  7.  Therefore,  referring  to  equations  (15),  and  letting  circled 
numbers  refer  to  columns  in  table  I,  we  have  that 


q.li  =  Z  ^3)  =  0.168 


^•12  “  0-21  ~  ^3^  X 

©  -  0.006 

^13  =  <131  *  Z  ©  X 

@  =  0.228 

*14  =  *41  “  -S  ©  >< 

@  =  -0.091 

<12 2  =  2  ©2  «  0.212 

0-23  =  1g2  =  -Z  ©  X 

©  =  0.117 

9-24  ~  ^42  *  Z  ©  X 

©  =  0.389 

^33  =  2©2=  0.415 

0-34  =  14g  =  -Z  ©  X 

©  =  0.051 

^44  =  =  °*985 

Inserting  these  numbers  into  the  expression  (l8)  for  D,  we  obtain 


Furthermore, 


0.168 

0.006 

0.228 

-0.091 

.006 

.212 

.117 

.389 

.228 

.  .117 

.415 

.051 

-.091 

.389 

.051 

'  .985 

0.000343 


0.212 

0.117 

0.389 

.117 

.415 

.051 

=  0.0145 

.389 

.051 

.985 

0.168 

0.228 

-0.091 

.228 

.415 

.051 

=  0.0115 

-.091 

.051 

.985 

3E 
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’m 

0.168 

0.006 

-0.091 

d3  = 

.006 

.212 

.389  =  0.00744 

* 

•• 

-.091 

.389 

.98  5 

0.168 

0.006 

0.228 

d4  - 

.006 

.212 

.117  «  0.00177 

.228 

.117 

.415 

Substituting  these  numbers  and  the  value  of  M^0)  given  in  equations  (25) 
into  equation  (20),  we  obtain 


|  AZ  |=  0.194 

|az'|=  .173 

|AP  |=  .139 
j  A  P  *  j  =  .068 

where  A 2  denotes  the  error  in  z,  A2'  in  i',  etc.  On  a  percentage 
basis,  this  implies  that  2,  2',  andp',  respectively,  are  known  to 
within  about  l4,  6,  23,  and  33  percent. 

Calculation  of  the  error  in  the  stability  parameters. -  It  has  been 
shown  how  the  errors  in  the  parameters  1,  1' ,  3,  and  3*  may  be  found. 
However,  the  problem  of  interest  to  the  aerodynamicist  is  the  calculation 
of  the  errors  in  the  stability  parameters  b,  k,  Cx,  and  CD  of  equa¬ 
tion  (l).  While  q(t)  can  be  written  directly  as  a  function  of  b,  k,  Ci, 
and  C0  instead  of  the  parameters  2,  2’,  3,  and  3’  and  the  method 
described  above  applied,  it  is  believed  that  the  following  general  con¬ 
siderations  simplify  the  calculations. 

Suppose,  as  before,  that  q(t)  is  a  function  of  the  m  param¬ 
eters  xjj,  k=l,...,m,  and  that  the  errors  Axjj  in  the  xj^  have  been 
found.  Suppose  further  that  there  are  n  other  parameters  y,,...,y  y 
each  of  which  is  a  function  of  xx,  x2, ...,  and  xm.  Thus,  11 


9 


'N 


y j  =  yj(xij***j xm) y  j - 1 , • • • , n 

The  error  Ayj  in  yj  may  then  be  estimated  from  the  formula 


It  should  be  noted  that  equation  (28)  gives  a  quite  pessimistic  value 
of  the  errors. 
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From  the  definition  of  2  and  2  ’ ,  we  may  now  write 


b  = 
k  = 


-22 

22+2* 


and  so. 


Also, 


8b 2  8b  _  8b  _  8b 

82  ,82»  8p  83* 


—  =  *2 
82 

—  =  22' 
82' 


M-  =  0 
83*  . 


(2  9) 


Thus,  applying  equation  (28)  using  the  error  values  (27),  we  see  that 

|Ab  |  =  2  |  A  2  |  =  0. 388 

|Ak  |  =  2  |  2A2  |+2  |  z»A2»  |  =  1.59 


and  since,  from  equations  (29), 

b  =  2.732  \ 
k  =  11.30  J 

it  follows  that  b  and  k  are  known  to  within  about  l4  percent. 

This  same  scheme  will  now  be  applied  to  the  computation  of  the 
errors  in  Cx  and  CQ.  It  will  be  noticed  that  since  Ci  and  CD  are 
quite  complicated  functions  of  2  and  2',  our  task  will  be  considerably 
more  difficult  than  it  was  when  computing  |Ab|  and  |Ak|. 

As  before,  suppose  the  input  S(t)  is  zero  for  all  t  greater 
than  ta.  Make  the  following  definitions: 
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/T-'Ct 

e"^T  5(t)  cos  2’ TdT 
*i 

/^a 

e"^T6(T)  sin  2 * tcLt 
ti 
ta 

Ri  =  /  Te"^T6(T)  cos  2’TdT 

ra 

Ri*  =  /  Te"^TB(T)  sin  2'tcLt 


so  that 


3Si  SSi' 


-  -  Ri 


Also  define 


dSj  dSi*  ^  t 

d2»  =  hi  =  _  Ri 


Oq^  =  e^^(Si’  cos  2'ti  -  Si  sin  2’ti) 

aii  =  -eZti[(2,Si-2Si')  cos  2'ti  +  (2Si+2'SiT)  sin  2’ti] 

ai  =  2' [  q.e(ti)  -  e2ti(p  cos  2’ti  -  p*  sin  2'ti)] 


IN 

■I 


n 

-I 


voi^il  ao  =  )  ctoiai 


a 

-Z 


aio  -  7  ^i0©!  an 


-Z* 


aliai 


With  these  definitions,  it  is  shovn  in  reference  2  that  Co  and  Cx 
satisfy  the  equations 


aoo  C0  +  a0j_  C-l  =  &o 
aio  Cq  +  an  Cj_  =  ai 
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Therefore,  if  P  denotes  any  one  of  the  four  parameters  l,  l1,  3>  3*> 
the  derivatives  8c0/Sp  and  can  he  found  from  the  equations 


8c, 

_  SaQ 

_  8a0o 

5a0i 

8P  +  &01  8P 

8p 

C°  A 

op 

1  ^P 

SCi 

Saj 

8a10 

^an 

8P  +  &11  8P 

8p 

C°  8p 

Cl  ^P 

The  problem  before  us  is  the  calculation  of  these  derivatives. 

From  the  definitions  (32)  and  from  equations  (3l),  it  follovs  that 
=  e2ti[(tiSi»  -  Ri’)  cos  Z'ti  -  (tiSi  -  Ri)  sin  Vt±) 

8  z 

—  =  tjttii  +  Ooi  +  eUi[(z«Ri  -  ZRi’)  cos  l't±  +  (ZRi  +  Z’Ri')  sin  Z'ti] 
8z 

5^1  =  -  Z  *t-,elti(p  cos  Z'ti  -  p*  sin  l't±) 

8  z 

=  eZtl  [(Rq  -  tqSi)  cos  Z’tq  +  (Ri‘  -  tiSq')  sin  Z’tJ 
8Z' 

^  =  -eZti{[Si  +  z(tiSi  -  Ri)  +  Z'(tiSi'  -  Ri’)I  cos  Z’ti  + 

8z ’  t  \ 

[Si*  -  Z'(tiSi  -  Ri)  +  Z(tiSi'  -  Ri’)]  sin  Z'ti 

— ^  =  —  +  z»tieZti(p  sin  Z'ti  +  3'  cos  Z'ti) 

8z«  V  y 

^oi  =  Q 

83 

!iii  =  o  > 

83 

ccs  !'ti 

a?  J 


^ai  _  „  .  Zti  .  , ,  j. 

- -  =  z'e  1  sm  Z'ti 
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The  quantities  e~lt5(t)  cos  Vt,  e“ltS(t)  sin  Vt,  te_it6(t)  cosZ’t, 
te"^5(t)  sin  Z’t  must  "be  integrated  to  find  S^,  Si’,  Rj^  and  Rq* . 

These  quantities  are  calculated  in  table  II  and  are  plotted  versus  time 
in  figure  1.  They  were  integrated  by  means  of  a  planimeter  to  obtain 
the  values  displayed  in  table  II.  Referring  to  table  II  and  equa¬ 
tions  (33)>  it  may  be  seen  that 


a00  =  2  @2  =  0.000399 

aoi  -  a1Q  =  Z  @  X  © 

an  =  r  ©  =  0.0241 


-0.00292 


a0  =  2  @  X 
ax  =2  ©  X 


-0.0615 

0.494 


and  that 


8a, 


00 


8z 


=  2Z 


X 


da01  =  8a1Q  _ 

8z  bl 

ba  11 


bl 

dap 

bl 

8ax 

Hi 


=  22 


=  2 


X 


=  2  GO)  X 


=  -0.000152 


+  2  GO)  x 


=  -0.01264 


+  ZGb)  x 


+  266)  X 


=  0.001388 


*  0.00772 


-  -0.0851 


8a, 


■00 


bl' 

8apj_  8a10 


bl' 

^an 
bl' 
8a, 


22(23)  X  ^9)  =  0.000232 
2^  X  +2 


8  2* 
=  22 


-  0.01226 


=  -0.001571 


^  =2@  x  ©  +2  ©  x  ©  =  -0.0381 
~  =  2  ©  x  ©  +  2  ©  x  © 


=  0.280 
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The  large  errors  in  some  of  the  parameters  are  worthy  of 
special  comment.  It  should  he  noted  first  of  all  that  the  particularly 
large  error  in  CD  is  completely  consistent  with  what  has  been  found 
for  some  time  by  empirical  means.  It  has  been  found  that  upon  repeat¬ 
ing  a  set  of  flight  records,  making  the  two  sets  identical  as  far  as 
that  is  possible,  two  entirely  different  values  of  CQ  have  often  been 
obtained. 

Secondly,  it  will  be  recalled  that  in  the  derivation  of  the  formula 
for  the  errors  it  was  required  that  the  errors  be  small;  this  is  cer¬ 
tainly  not  the  case  for  some  of  the  errors.  For  this  reason,  the  error 
of  37-5  in  C0  particularly  may  not  be  considered  as  definitive.  The 
theory  developed  herein  fails  for  such  a  large  error,  and  the  actual 
error  in  CQ  may  be  much  larger  than  that  calculated.  However,  this  is 
not  serious,  for  an  error  as  large  as  99  percent  renders  the  calculated 
parameter  value  meaningless  in  any  case,  and  it  matters  not  at  all 
whether  the  error  is  99  or  199  percent. 

Finally  and  most  important  is  the  following  conclusion  regarding 
the  entire  philosophy  of  calculation  of  stability  parameters  from  pitch¬ 
ing  velocity  data  alone  in  response  to  an  elevator  pulse:  This  experi¬ 
ment  is  ill  designed  for  the  calculation3  of  Co.  In  general,  when  a 
pulse  is  applied  to  produce  a  set  of  pitching  velocity  data,  no  more 
should  really  be  expected  from  the  analysis  than  the  period  and  the 
damping  parameters  k  and  b.  If,  in  addition,  the  analysis  results  in 
a  value  of  C-^  whose  error  is  within  reasonable  limits,  this  should  be 
regarded  as  fortuitous.  It  should  be  stated  that  no  example  has  yet 
been  found  for  which  CD  may  be  calculated  at  all  accurately. 

Such  a  negative  comment  as  the  preceding  deserves  a  remark  on  the 
possibility  of  the  most  accurate  calculation  of  Co-  If  one  were  to 
measure  a  step  response  rather  than  a  pulse  response,  the  analysis  of 
these  data  would  surely  lead  to  more  reliable  values  of  C0.  The 
reason  for  this  lies  in  the  easily  proved  proportionality  of  the  steady- 
state  value  of  q(t)  with  CD. 


Q 

This  conclusion  is  drawn  for  the  parameter  CG  alone  and  not  also 
for  Cj_,  as  would  seem  to  be  indicated  by  the  error  of  72  percent 
in  Clt  primarily  because  the  error  in  Ci  is  not  always  as  large  as 
this.  With  some  data  which  have  been  analyzed,  the  error  in  has 

been  far  less,  o'f  the  order  of  30  percent. 
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CONCLUDING  REMARKS 


A  formula  las  been  given  which  may  he  used  to  find  the  errors  in 
the  parameters  obtained  from  a  curve -fitting  process.  The  method  of 
derivation  did  not  use  the  concepts  of  probability  theory,  since  the 
latter  ideas  lead  to  quantities  for  which  there  is  no  known  method  of 
solution.  However,  the  result  that  the  formula  obtained  bears  a  close 
relationship  to  the  classical  probabilistic  formula  in  case  the  curve - 
fitting  problem  is  linear  (the  only  case  for  which  such  a  formula  has 
been  derived  and  proved  valid)  is  proved. 

As  may  be  seen  by  studying  the  example  given,  the  ease  of  appli¬ 
cation  of  the  method  is  directly  dependent  on  the  simplicity  of  the 
form  of  the  function  which  has  been  fitted.  Thus,  if  b,  k,  Cx,  and  C0 . 
are  the  aerodynamic  parameters  occurring  in  the  differential  equation  (l), 
the  errors  in  b  and  k  may  be  calculated  fairly  rapidly;  the  errors 
in  Cx  and  CG  require  more  time.  However,  if  the  method  of  least 
squares  (reference  2)  is  used  for  the  curve -fitting  problem,  a  great 
many  of  the  computations  which  are  needed  to  find  the  errors  will  have 
already  been  performed  in  the  process  of  finding  the  parameters. 

The  method  actually  weights  favorably  the  method  of  least  squares 
in  another  manner.  The  errors,  as  given  by  equation  (20),  are  propor¬ 
tional  to  the  square  root  of  the  sum  of  the  squares  of  the  residuals, 
which  quantity  is  minimized  by  the  least-squares  process.  However,  it 
is  not  believed  that  this  weighting  is  a  serious  limitation,  for  no 
matter  what  means  are  used  for  fitting  the  curve,  if  this  latter  method 
is  to  have  significance  and  is  to  lead  to  a  good  fit  of  the  data,  the 
sum  of  the  squares  of  the  residuals  will  also  be  near  its  minimum. 

Finally,  it  has  been  shown  that  pitching  velocity  data  in  response 
to  a  pulse  are  not  alone  adequate  to  compute  the  parameter  C0  (and,  to 
a  lesser  extent,  Cx)  occurring  in  the  differential  equation  of  motion. 
Usually,  all  that  may  be  obtained  from  such  data  are  the  parameters  b 
and  k,  which  determine  the  damping  and  the  period  of  the  oscillation. 

In  some  cases,  Cx  may  also  be  obtained  with  reasonable  accuracy. 

Ames  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Moffett  Field,  Calif.,  July  30,  1952 
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